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Abstract 

The bound states of the fermionic 3 He(2 3 Si)+ 3 He(2 3 P ; ) system, where j 



o 

(N 



0, 1,2, are investigated using the recently available 
ab initio short-range 1 ' J ' 3 £+ M and l - i - 3 Tlg M potentials computed by Deguilhem et al. (J. Phys. B: At. Mol. Opt. Phys., 2009, 
42, 015102). Single-channel and multichannel calculations have been undertaken in order to investigate the effects of Coriolis 
and non-adiabatic couplings. The possible experimental observability of the theoretical levels is assessed using criteria based 
upon the short-range character of each level and their coupling to metastable ground states. Purely long-range levels have been 
identified and 30 short-range levels near five asymptotes are suggested for experimental investigation. 
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1 Introduction 

Photoassociation (PA) of ultracold atoms, in which two in- 
teracting ultracold atoms are resonantly excited by a laser to 
bound states of the associated molecule, is a widely used tech- 
nique to study the dynamics of ultracold collisions in dilute 
quantum gases. Of particular interest is PA in metastable rare 
gases where novel experimental strategies based upon their 
large internal energy can be implemented. 

Photoassociation of ultracold bosonic metastable 4 He* 
atoms, 4 He(ls 2s 3 S), to excited rovibrational bound states 
that dissociate to the 4 He(ls 2s 3 S) + 4 He(ls 2p 3 P 7 ) limits, 
where j — 0, 1 , 2, has been observed by many groups. The ob- 
servations include over 40 states lying within 14 GHz of the 
j = 2 asymptote fl] E] [3], six states within 0.6 GHz of the 
j = 1 asymptote [3] and some purely long-range bound states 
within 1.43 GHz of the j = asymptote [4|. Theoretical anal- 
ysis of the 7 = long-range states using single-channel [5 1 and 
multichannel [6| calculations based upon long-range Born- 
Oppenheimer potentials constructed from retarded resonance 
dipole and dispersion interactions gave excellent agreement 
with the measured binding energies. Analysis of the other 
states had to await the availability of short-range ab initio 
1,3,5 E+ [( and 1,3 ' 5 n S)M molecular potentials JTHU and was ini- 
tially restricted to single-channel calculations [7 , 8 1 which ne- 
glect non-adiabatic and Coriolis couplings. Very recently a de- 
tailed theoretical analysis of the entire 4 He(ls 2s 3 S) + 4 He(ls 
2p 3 P y ) system has been completed [9j. The role of these cou- 
plings was investigated using single-channel and multichan- 
nel calculations with the input potentials constructed from the 
short-range ab initio potentials of Deguilhem et al. [ 8 1 matched 
onto long-range retarded resonance dipole and dispersion po- 
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tentials. The multichannel calculations also permitted criteria 
to be established for the assignment of the theoretical levels 
to experimental observations based upon the short-range spin 
character of each level and their couplings to the metastable 
ground states. Excellent agreement was obtained for the num- 
bers of observed levels and their binding energies after appli- 
cation of a 1% increase in the slope of the 5 £g„ and 5 n g „ 
potentials near their inner classical turning point. 

In contrast, PA of fermionic metastable 3 He* atoms, 3 He(ls 
2s 3 S), is relatively unexplored although they have been cooled 
and trapped [ 1 1 with comparable densities and temperatures 
to those of 4 He* atoms. The non-zero i = 1/2 nuclear spin 
of 3 He* gives rise to hyperfine structure with splittings com- 
parable to the fine structure splittings of 4 He* which has no 
nuclear spin. Consequently the patterns of energy levels is ex- 
pected to be quite different for the fermionic and bosonic sys- 
tems. A small number of long-range states in 3 He* has been 
predicted by Dickinson [ 1 1 1 but this was a single-channel cal- 
culation, thereby neglecting Coriolis and non-adiabatic cou- 
plings, using only long-range van der Waals and retarded res- 
onance dipole interactions. The availability of the short-range 
potentials of Deguilhem et a/.[8 1 now permits a detailed theo- 
retical investigation of the fermionic 3 He(ls 2s 3 S) + 3 He(ls 
2p 3 P ; ) system similar to that undertaken by Cocks et al.[9\ 
for the bosonic 4 He* system. 

In the absence of any observations of bound states in this 
excited 3 He* system, we present predictions as to which of 
our calculated bound states may be experimentally observable. 
We assume any experiment will use magnetic trapping of the 
3 He* atoms, requiring all atoms to be in the fully stretched 
low-field seeking / = 3/2, m/ = 3/2 magnetic substate of 
the metastable 2s 3 Si level in order to strongly suppress loss 
through Penning ionization. Consequently we assess the ex- 
perimental observability of each excited level in terms of its 
coupling to this state. In addition, we consider the likelihood 
of ionization losses from these excited levels due to inelastic 
collisions in the short-range region. 

Atomic units are used, with lengths in Bohr radii aq = 



1 



0.0529177209 nm and energies in Hartree Ef, 
27.21 1384 eV. 
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2 Theory 

2.1 Multichannel equations 

The formalism for the excited 3 He* system requires modifica- 
tion of that presented in Cocks et al.\9\ for the excited 4 He* 
system in order to include hyperfine structure. 

The total Hamiltonian for a system of two interacting atoms 
i = 1,2 with reduced mass fi, interatomic separation R and rel- 
ative angular momentum 1, which possess both fine structure 
and hyperfine structure is 



H T + H rot + H e i + Hf s + Hhf s 
where T is the kinetic energy operator 

h 2 d ( , d 

2fj,R 2 dR \ dR 
and H mt the rotational operator 

i 2 



H„ 



2\lR 2 ' 

The total electronic Hamiltonian is 

II,-. ll\ I "2 • 1 1 



12- 



(1) 



(2) 



(3) 



(4) 



where the unperturbed atoms have Hamiltonians Hi and their 
electrostatic interaction is specified by H\2- The terms Hf S and 
i/hfs in equation (Q~|i describe the fine structure and hyperfine 
structure respectively of the atoms. 

The multichannel equations describing the interacting 
atoms are obtained from the eigenvalue equation 



(5) 



for the total system by expanding the eigenvector in terms of 
an appropriate basis |<J> ( ,) = \<t> a (R,q)) where a denotes the set 
of approximate quantum numbers describing the electronic- 
rotational states of the molecule and q denotes the interatomic 
polar coordinates (9,(p) and electronic coordinates (ri,r2). 
Using the expansion 



i 



£-G fl (tf)|<t> Q ) 



(6) 



and forming the scalar product (<!>„/ |//| V F) yields the multi- 
channel equations 



E { T $a (*) + %'a (R) - E8 a , a ] G a (R)}=0, 



(7) 



where 



and 



V a >a {R) = («V I [^rot + H el + H fs + H Ms ] \<P a ) 



(8) 



(9) 



We assume the /^-dependence of the basis states is negligible 
so that the radial kinetic energy term is diagonalized: 



T c ? a (R) 



H 2 d 2 G a 
2fl dR 2 



(10) 



2.2 Basis states 



For two colliding atoms with orbital L,, spin S, and nuclear i, 
angular momenta, the unsymmetrized body-fixed states in the 
coupling scheme 

j ; =L,+S ; , fi=Ji + U, f = fi +f 2 , f = f + l (11) 
are (see appendix for details) 

\(Yih hh )a, {Yihhh)B,f, %, T, m T ) (12) 

where ji = {fi,Li,Sj}, fi representing any other relevant quan- 
tum numbers, and (A,B) labels the two nuclei. The projec- 
tions of an angular momentum J onto the space-fixed Oz and 
inter-molecular axis OZ with orientation (9, (p) relative to the 
space-fixed frame will be denoted mj and Q.j respectively. 

In order to construct states symmetrized with respect to the 
total parity P T we note that P T = P L P s PiX N where Pl,Ps,Pi 
are the inversion operators on the orbital, electronic spin and 
nuclear spin states associated respectively with 



L = Li+L2, S = Si+S2, 



= 11+12 



(13) 



and Xn permutes the nuclei labels. The states of total parity 
are then (see appendix) 

\{a 1 )A > (a 2 )B,f^,T,m T ;PT) = 
Npt [\{ai) A ,(a 2 ) B ,f,<l>J,mT) 
+ PTP l P2(-l) f - T \{a l ) A ,{a 2 )B,f,-4>J,m T )] (14) 

where a,- = {ji, ji , ii , fi }, Pi = ( — 1 ) L ' 1 is the parity of the atomic 
state and</> = |£2/| = |£lr|- The normalization constant 

is Np T = 1 /^2(1 + 5^,0). For = equation (TBI i gives the 
selection rule PrAft(-l) /_r = 1- 

The states symmetrized with respect to Xn are (see ap- 
pendix) 



\a u a 2 ,f 1 <j).T,m T ;P T ,XN) = 

Nx N [\{ai) A .Xa 2 ) B JA-T 1 m T \P T ) 
+ £n I {Oq)a , («i )b , /, <P , T, m T ;P T ) 



(15) 



2 



where Pn = (— l) 2 ' 1 indicates bosonic or fermionic nuclei 
(where i\ = i 2 is assumed), Nj is the number of electrons on 
atom i, the normalization constant Nx N is 1 / y/2(l + 8 ai .a 2 ) 
and the phase factor is 



e N = P N P T PiP 2 (-l 



Jl+f2-f+NiN 2 



(16) 



For (Xi = a 2 , equation (T5[ gives the selection rule 

p N p T PiP 2 {-iy i+h - f+NiN2 = 1. 

It is convenient to introduce the simplified notation 



\a u a 2 ,<p) = Nx N [\(ai)A,(a 2 )Bj,<t>,T,mT) 
+e N \{a 2 ) A ,(ai) B J,(j) 1 T,m T )] 

so that the states ( TlSb can then be written 



(17) 



\ai,a 2 ,f,<j),T,m T ;P T ,X N ) =7% [|oci,a 2 ,0) 

+ P T P l P 2 (-iy- T \a 2 ,a l ,-$)]. (18) 

The eigenstates of H e \ are the body-fixed states arising from 
the couplings L = Li + L2, S = Si +S2 and must be symmet- 
ric under the action of PlPs- 

\Yxyi,LSQ. L Q. s ;w) =N w \\{Y\)a{Yi)b,LSQ.lQ. S ) 

+e w \{yi)A(Y\)B,LSa L Q. s )] (19) 



where N w = 1/^/2(1 ■+Sy^), w = 0(1) for gerade (unger- 
ade) symmetry and 



e. = (-i) 



W+Lj +L 2 +Si +S 2 - S+N l N 2 



PlP 2 - 



(20) 



The relationship between the two bases (TTTb and ( fT9l is ob- 
tained using (see appendix) 

\a u a 2 ,(j)) = 

N XN N H .\T,m T ,t) £ £ ^f^^S 
ijLia /.su,n„ 

x [(\g) + \u)) + e(\g) - \u))]\(h) A (h) B ,i£li) (21) 

where the coupling coefficients Fj^ 2 ^ and F^q^ are given 
in the appendix (the quantum numbers (L,-,S,-, !,•) have been 
suppressed) and we have introduced the notation 

\g) = \jxy 2 , LSQ. L Q. s ;g) , \u) = \ nn, LS£l l £l s ; u) (22) 

for the eigenstates of gerade and ungerade symmetry. The 
rotational states are 



\T,m T ,4>) 



IT + 1 

-^—D T m ; A (<p,e,o), (23) 



where Df n * ^(^,0,0) is the Wigner rotation matrix, and the 
phase factor is 



e = P N P T (-l) 2i ^ +i - 2s . 



For the 3 He(ls 2s 3 S) + 3 He(ls 2p 3 P ; ) system, <X\ = 
(7 1 ,0,l,l,l/2,/ 1 )anda 2 = (72,l,l,/2,l/2,/2) and (ED re- 
duces to 

\a u a 2 ,<p) = 

\T,mj, <p) £ £ (-l) l - J2 [UM2Sj 2 } l/2 



ijif r \Sj 



XC^C^\ 1/2 1/2 \){\)s 

x l -[(\g) + \u))+p T (-iy(\g)-\u))] 

x\(h)A(h)B,i&i) 



where Cm/mim is a Clebsch-Gordan coefficient, 



(25) 

a b c 
d e f 



1 "'2 
a b c 

and ^ d e f ^ sre Wigner 6 — j and 9 — j symbols respec 

g h i 

lively and [ab . . .} = (2a + 1) X (2b + 1) x .... 
2.3 Matrix elements 

The multichannel equations (O require the matrix elements of 
H TO t, H e \, Hf s and Hhfs in the basis ( TTTl ). Using the notation 
I a) = |3> a where a = {a u a 2 J,(j>,T,m T ,P T ,X N } then 

the rotation terms are 

(a'\l 2 \a) = h 2 8 p , p ,{[T(T + l)+f(f+l)-2<j> 2 ]8 r ^ 

~~ K Tf<j) ~ K Tf<j) "V^+l} ' (26) 

where the Coriolis coupling terms are 



K. 



Tf<P 



[r(r+i)-M±i)]s 

x[/(/+l)-0(</>±l)]? 



(27) 



(24) 



and p denotes the set of quantum numbers 
{a u a 2 ,f,T,m T ,P T }. 

The electronic matrix elements can be expressed in terms of 
the Born-Oppenheimer (BO) molecular potentials 2s+l A°(R), 
where A = \Q.l\ and a is the symmetry of the electronic wave 
function with respect to reflection through a plane containing 
the internuclear axis, using 

H el \ 7l72 ,LS£L L £L s ;w) = f^A^tf) \YiY2,LS£l l £l s ;w) 

(28) 

where is the asymptotic energy of the state. The result is 
(see appendix) 

(a'\H el \a) = 

Kn> I I (-l) J2+ ^ +; ' +/ [/i/i72/^2; 2 //'] 1/2 
jfisa L o.i 
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+P T {-iy [ 2s+l A+(R)- 2s+l A+{R)] } (29) 

where ( ° * £ J is a Wigner 3 — jf coefficient, £2$ = 

(j> — Q.i — Ql and tj denotes the set of quantum numbers 
{ 7i j 72 > , 017- ,Pt}- This equation differs from that given by 
Dickinson ifTTI by an overall phase factor (— 1) 1_!_£2 ' and the 
phase of the A g — A u term. 

The matrix elements of the fine structure and hyperfine 
structure are best expressed in the basis 

rrij-mi. m^.m^. 

(30) 

For convenience we omit the label nif i from these states as the 
matrix elements of Hf s and //hfs are independent of due 
to rotational invariance. We assume that the fine structure is 
independent of R and exclude couplings to the singlet atomic 
state S: = so that its contribution is 



K|# fs |ot;) = 8^, AE { *j.. 



(31) 



The fine structure splitting AEj^ for the 2s 3 Si level vanishes 
and the splittings AE* ■ for the 2p 3 Po and 2p 3 Pi states rel- 
ative to the 2p 3 P 2 level are 31.9088 GHz and 2.2922 GHz 
respectively [ 13]. 

Matrix elements for the hyperfine structure have been ob- 
tained by Hinds et al.\ 12] and Wu and Drake [13]. We choose 
to use the expression of Wu and Drake but exclude couplings 
to the Si = atomic states. The matrix elements are therefore 



(a-|ff h f s |0!;; 



W. 



Ufi 



C Si V6(-l)W>X Si 



_2) Sj (_i)J»+S/+Jf 



°i Ji 

ji Si 



Li fi S, 
ji ^ 1 



Li 1 U 
-MOM 



-1 



Li Li 2 
Si S, 1 



Ji Ji 



1 



Li 2 Li 
-MOM 



(32) 



where these expressions are to be evaluated with M = Li, 



ji U J V 



-1 



(33) 



and 



^ = -(25 i+ l){ f ^ f}. (34) 
The hyperfine structure parameters (in MHz) are [13] 

Ci = -4283.85, D\ = -28.145, £1 = 7.126. (35) 



The inclusion of hyperfine structure using (132t couples states 
with the same L,, Sj and /, but different ji and for the He 2p 3 P 
manifold the states (j,f) — (0,1/2) and (1,1/2) are signifi- 
cantly coupled as are the pair (j,f) — (1 , 3/2) and (2,3/2). The 
eigenvalues of fff S + i/hfs give the following energies for the 
hyperfine levels expressed relative to the state j = 2,f = 5/2: 
0, 1780.851, 6292.906, 6961.065 and 34385.941 MHz. The 
eigenvectors give the mixing coefficients which are then used 
to modify the purely algebraic transformation to the hyperfine 
case given by equation (13 Oi l. For the hyperfine structure of 
the 2s 3 S level we adopt the splitting of 6739.701177 MHz 
as measured by Zhao et fl/. lfl41 . This data then gives the 
ten asymptotic energies of the separated pairs of atoms 
as 0, 1780.851, 6292.906, 6739.701, 6961.065, 8520.552, 
13032.607, 13700.766, 34385.941 and 41 125.642 MHz. 

We assume that the fine- and hyperfine- structure of the in- 
dividual atoms is not affected by their participation within the 
dimer, so that we may write 

(a'\Hf s +H Ms \a) = 

8 a , a (AE* h +AEl h ) + 8 a , a (« |#hf S |«i) 

8cd,,a,2 

+(a^ hfs |a 2 )5 a ;, ai ) (36) 

where a denotes the set of quantum numbers 
{f,(j),T,m T ,P T }. 

The total matrix element V a i a (R) is therefore diagonal in 
{T,Pj} and independent of mj. The mj -degenerate discrete 
multichannel eigenenergies of (|7]i are then £r,/V;v where v la- 
bels the rovibrational levels. 

2.4 Single-channel approximation 

The single-channel approximation involves the neglect of the 
Coriolis couplings in (|26| | and non-adiabatic couplings in the 
kinetic energy term. At each value of R the single-channel 
potential is formed by diagonalizing the matrix: 



Vi = if!\H A \a) + (a'|(# fs +# Ms )|a> 



a'\?\a) 
2jxR 2 



(37) 
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where (a'\P\d)s is the part of (f26t diagonal in <j). The corre- 
sponding ^-dependent eigenvectors are 



Y^C an {R) 



(38) 



and the adiabatic potential is given by Vf\R) = 
Y^a'aC^ly^ijCan- Each channel \n) can be labelled with 
the notation { <p , T, mj ,Pt}- 

The adiabatic eigenvalue equation for the rovibrational 
eigenstates |y/„. v ) = R~ l G„^,(R)\i), where n = {(j>,T,niT,PT}, 
is then obtained by neglecting the off-diagonal (non-adiabatic) 
couplings between different single-channel states in the ki- 
netic energy term so that 



{n'\ijG n , v {R)\n) 



h 2 d 2 G,,, v 
2jiR dR 2 



(39) 



The radial eigenvalue equation for the rovibrational states is 
then 



h 2 d 2 



■ E n-V 



G„, v (R) = 0. 



2.5 Input potentials 

The required Born-Oppenheimer potentials > 3,5 £;j 



(40) 



and 



1,3,5 



Tlg. u were constructed as in Cocks et a/.|9) by matching 
the ab initio short-range potentials of Deguilhem et al. [8 1 onto 
the long-range dipole-dipole plus dispersion potentials 



V l A ong (R) = -f 3A (R/k)C 3A /R 3 - C 6A /R 6 



CfjR*-C 9A /R»-C W A/R 10 , 



(41) 



where is an R- and A-dependent retardation correc- 
tion lfl6l . k = X/(2n) = 3258.12oo where A is the wave- 
length for the 2s 3 S-2p 3 P transition and the parameters C„a 
were taken from Zhang et a/. lfl5ll . 

Motivated by our study of the 4 He* system @, we choose 
to vary the quintet potentials through a modification of the 
slope of the potential at the inner classical turning point by 
introducing a multiplicative factor c through the smoothing 
function 

( V{R){\+2c) R<Ri 
V'(R) = < V(R)[l+c(l+cosa(R-Rx))] Ri<R<R 2 , 
{ V{R) R>R 2 

(42) 

where R\ = 5ao, R 2 = 10 ao and a = nj (R 2 — R\). The value 
c = 0.005 represents a 1% variation which is quickly turned 
on through the region 5 < R < 10ao. Its effect is to deepen 
the minimum of the attractive 5 Yl g potential at R = 5.387 ao 
by 0.985% and move it to a smaller interatomic separation by 
0.003 ao- The depth of the minimum in the 5 £+ potential at 
R = 6.268 ao is increased by 0.851% and is moved towards a 



smaller separation by 0.010 ciq. The other quintet potentials 
5 S+ and 5 n [( are not significantly affected as they are repul- 
sive. 



3 Results 
3.1 Calculations 

The coupled-channel equations ^ and the single-channel 
equation (l40l are of the form 



G(R)=0, 



(43) 



where for the case of coupled-channels, G is the matrix of 
solutions with the second subscript labelling the linearly in- 
dependent solutions. These equations were solved using the 
renormalized Numerov method ifTTI with the eigenvalues of 
the purely bound states determined by counting the nodes 
of the determinant |G(/?)| and the energies of resonances ly- 
ing within open channels by using a search procedure based 
on Cauchy's argument principle applied to the determinant 
D{E) = |R,„ — R m +il where R„, and R,„+i are ratio matri- 
ces for the outward and inward integrations respectively of the 
renormalized Numerov method. Further numerical details are 
given in Cocks et q/.||9l. 

3.2 Observability Criteria 

In order to predict the likelihood that calculated bound levels 
may appear in future experiments, several properties are de- 
termined for each bound level or resonance that we isolate. 
The simplest of these is the proportion P s h rt of wave function 
present at close interatomic distances, defined as R < 20 ao 
and henceforth referred to as the short-range region. This 
property is extremely useful in classifying results since ion- 
ization losses, which arise from the inelastic collisions 



He*+He* 
He*+He* 



He + He + +e 
Het + e~ 



(44) 



only occur in the short-range region. As has been observed 
in bosonic metastable helium, there exist indications of purely 
long-range states in the fermionic dimers investigated here, 
and we define these by P s hort < 10~ 10 . 

If the level extends into the short-range region then an in- 
dication of its propensity for ionization is obtained from the 
proportion P str of wave function that is in the spin-stretched 
S = 2, ; = 1 configuration: 



ab 



^La.bPab 



(45) 
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where 



20 Oq 



P ab = (a\b) / G a (R) dR 



(46) 



and \b) = \ Ji Yz, LSQ.lQ.sw) \ (z'i )a (h)B> iQ) is tne complete LS 
basis state. The transformation between the bases used here 
can be found from equation (fSTJ. As in the 4 He* case, the 
the ionization rate of the dimers is significantly reduced in the 
spin-stretched state fl8l . Hence, a large proportion of wave 
function in the spin-stretched state is essential for the level to 
have a lifetime long enough to be observed in experiment. 

Finally, for a resonance to be observed in PA experiments, 
it must be strongly coupled by a laser pulse to the metastable 
manifold 3 He(ls 2s 3 Si)+ 3 He(ls 2s 3 Si). For radiation of cir- 
cular polarization ex the coupling between a metastable dimer 
state and the excited dimer state is due to the interaction 
Hm ~ £a ■ d where d is the molecular dipole moment and is 
given by 

(e'\H- mt \g) = 



■ ( ,a 1 / 2T + 1 TlT , T1T , Nx N 

[ ' Vi^V2r 7 TT Cm ^^ L w~ al 

xi l (-iy-His]i/f[fy' 2 m l/2 

iSj' il.'AAl' 

Yr fif r W r Sif / 1 1 / 2 1 

xl -a / a i a / , u /3a s n / L 'n5a,(|) | \ f s ] 

1 k /If 1 15 

x { 1/2 1/2 i } { 1/2 1/2 i } (47) 

fx A 



A A f J I h h f 



where \g) 
ing to 



and \e) are basis states correspond- 
the 3 He(ls2s 3 Si)+ 3 He(ls2s 3 Si) and 
3 He(ls 2s 3 Si)+ 3 He(ls 2s 3 P 7 ) manifolds respectively, 
/5 = <p' — and only those matrix elements with P' T Pr = — 1 
and |/5 1 < 1 are non-zero. The atomic dipole moment is 
given by d at . Note that the metastable spin-stretched state 
has Pt = +1 symmetry and can therefore only be coupled to 
excited dimers of Pt = — 1 symmetry. 
Finally, we define the quantities 



a\H- m \g') / Gg,(R)G a , v (R) dR, (48) 



where g' enumerates all of the N„i spin-stretched metastable 
dimer states, and 



^full 



1 



g ga 



Y,(a\H m t\g) G s {R)G aiV {R)dR (49) 



where g enumerates all of the N g metastable dimer states. The 
true metastable radial wave functions G g (R) depend upon tem- 
perature, but in order to extract a single parameter for the ob- 
servability criteria, we take G g {R) — 1 as was done in the 4 He* 



case. This is valid up to a constant factor when the metastable 
wave functions do not change significantly with temperature. 
Although we focus on predicting resonances observable from 
experiments prepared with spin-stretched states in this pa- 
per, due to the overwhelming benefits from reduced trap loss, 
whenever it is convenient we also include the likelihood for 
couplings from other metastable states. Spin-stretched exper- 
iments are best described by the criterion g/ s ti, whereas ex- 
periments that do not polarize the metastable gas are best de- 
scribed by the criterion ^iiii- 

3.3 Single-channel 

The binding energies of long-range states obtained using a 
single-channel calculation are listed in Tables [T] and [2] The 
single-channel levels are labelled in terms of {T, <j>^}. 

Levels which are strongly coupled to the spin-stretched 
metastable dimer states are indicated by a superscript 1 . In the 
absence of existing experimental data, we use the same crite- 
rion to that obtained for the 4 He* case, that is £/ stI > 0.9 £/,. 
As these are long-range levels, there is no possibility of ioniza- 
tion and we can ignore the P str condition. Furthermore, levels 
that are strongly coupled to the unpolarized metastable dimer 
states are indicated by a superscript 2, where the criterion is 
s/ m > 0.9 E h . 

Of the 159 long-range levels found, 15 have a strong spin- 
stretched coupling, and 69 have a strong unpolarized cou- 
pling. In addition, there are 151 levels that possess some 
short-range character, and also satisfy the observability cri- 
teria. Some of these are very strongly coupled to the spin- 
stretched metastable state. However, we do not observe these 
levels once non-adiabatic and Coriolis couplings are turned on 
and so conclude that these levels are unlikely to be observed 
in experiment. 

3.4 Multichannel 

With all couplings included in the calculation, only those lev- 
els beneath the lowest asymptote are true bound states. In 
contrast to the situation in 4 He*, most of the levels lie above 
the lowest asymptote and, due to couplings to open channels, 
these higher lying levels almost always acquire a finite lifetime 
due to predissociation. These resonances possess complex en- 
ergies, where the imaginary component represents the reso- 
nance width, and are more difficult to isolate. As our search 
routine based on Cauchy's argument principle requires many 
solutions of the differential equations (0, we restrict the pre- 
dissociation width to be less than 100 MHz and only search 
within 2 GHz of the asymptotic energies that result from diag- 
onalization of the hyperfine structure. Additionally, we match 
only at two points, 100 and 300 ao, which may exclude a few 
levels from our search, although it can be argued on the ba- 
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Table 1 Single-channel rovibrational binding energies, in units of MHz, of long-range 0* and 1* states in 3 He(2 3 Si)+ 3 He(2 3 Pj). Energies 
given are relative to the energy of the specified asymptote. The superscripts 1 and 2 indicate those states which satisfy the strong coupling 
conditions s/ sts > 0.9 £/, and i<4ull > 0.9 E/, respectively. 



Symmetry 


State No. 


Asymp. No. 


v/T 





1 


2 


3 


4 


0+ 


5 


3 





904.113 




823.629 




639.703 








1 


183.326 




127.625 




9.53314 




6 


4 





347.642 




262.537 2 




75.7138 2 








1 


10.8473 












10 


6 







1422.19 




1278.29 










1 




467.747 




366.810 










2 




62.6480 




11.3740 






ll 


8 


16 




52.8477 2 




27.0296 










17 




11.6235 










12 


9 







202.645 2 




52.0862 2 










1 




13.4273 








(T 


7 


5 





6.10337 












10 


7 





374.065 


1425.59 


296.578 2 


1271.97 2 


126. 120 2 








1 


40.4103 


579.160 


8.45110 


449.478 










2 




173.330 




99.6096 










3 




37.9230 










11 


8 


16 




940.3 15 2 




815.427 2 










17 




503. 466 2 




417.433 










18 




227.478 2 




172.257 










19 




86.2920 2 




55.4097 










20 




24.2101 2 




9.30074 






12 


9 


16 




500.709 2 




319.184 2 










17 




131.414 




47.6889 










18 




19.7607 










13 


10 







741.860 2 




547.258 2 










1 




233.251 




130.685 2 










2 




51.1549 




11.7006 










3 




5.72522 








-i -I- 

l+ 


17 


7 







1269.28 


1214. 50 2 


1132.90 


1025.19 








1 




405.296 


366.915 


310.695 2 


238.337 








2 




74.8580 


55.8199 


29.9220 


1.87478 








3 




6.48198 










18 


8 







918.791 2 


869.335 2 


796.261 2 


701.007 








i 
i 




AH fxll'Z 
HOZ.O / / 


A0 1917 


J^J. 100 










2 




180.375 2 


160.938 


133.043 2 


98.3456 








3 




62.5650 2 


52.0562 


37.1800 2 


20.8902 








4 




15.6929 


10.8815 








20 


9 







438.234 2 


366. no 2 


264.239 2 


141. 524 2 








1 




106.356 


73.1955 2 


31.0864 










2 




14.2935 


5.04058 








22 


10 


5 




91.4710 2 


58.1199 2 


16.2073 2 










6 




8.96832 


1.27172 






l- 


11 


4 







526.589 1 


465.05 1 1 ' 2 


374.249 1 ' 2 


256.239 








1 




31.3952 


6.28362 1 * 2 








12 


5 







30.9200 1 ' 2 










17 


7 







342.560 2 


290.73 1 1 - 2 


214.234 1 ' 2 


114.816 2 
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Table 2 Single-channel rovibrational binding energies, in units of MHz, of long-range 2 and 3* states in 3 He(2 3 Si)+ 3 He(2 3 Pj). Energies 
given are relative to the energy of the specified asymptote. The superscripts 1 and 2 indicate those states which satisfy the strong coupling 
conditions ,s/ str > 0.9 Ef, and &?m\ > 0-9 E/, respectively. 



Symmetry 


State No. 


Asympt. No. 


v/T 





1 


2 


3 


4 


2+ 


6 


3 









1263.59 


1167.37 


1041.69 








1 






483.739 


400.323 2 


292.599 2 








2 






24.8937 








7 


4 









1119.79 2 


1013.79 2 


875.239 2 








1 






436.459 2 


376.891 2 


301. 106 2 








2 






143.028 


105.268 


58.7279 2 








3 






10.5484 








8 


4 









524.832 2 


441.585 2 


332.271 2 


2" 


12 


8 









932.947 2 


853.811 12 


749.864 1 ' 2 








1 






325.882 2 


278.283 2 


217.474 1 - 2 








2 






95.8030 


71.0373 


41.3769 1 - 2 








3 






19.9633 


10.0466 






13 


9 









87.8731 1 ' 2 


35.2247 1 ' 2 










1 






5.25868 






3+ 


4 


4 











1623.55 


1470.74 








1 








568.235 


468.826 








2 








137.420 


81.9036 








3 








0.85611 




3" 


4 


4 


16 








643.992 12 


551.467 1 ' 2 








17 








150.1 16 1 


97.8055 








18 








1.20481 





sis of spin-conservation of the laser coupling that resonances 
which exist solely inside this distance will very likely ionize 
and hence will not be observed in experiment. 

Beneath the lowest asymptote we find bound levels with 
only very weak coupling strengths. We therefore focus on the 
resonances that were successfully isolated. As these levels 
are not purely long-range, we must also consider the effect 
of ionization which reduces the level's lifetime and hence ob- 
servability. In our previous investigation of 4 He* we imposed 
a criterion of P str > 87.5%. However, although a large num- 
ber of resonances were found in 3 He* using the above method, 
very few satisfy the same observability criteria as 4 He* . In Ta- 
ble[5]we instead list the 30 resonances that are most likely to be 
observed in experiment, grouped by the nearest fine-structure 
asymptote. 

In contrast to the purely long-range levels in the 0+, J = I 
potentials of 4 He*, we do not find any single-channel long- 
range bound levels in the 3 He* potentials that remain bound 
after the inclusion of couplings to all accessible states, nor do 
we find any multichannel levels that can be described purely in 
terms of single-channel potentials. Again we must emphasize 
that the relative coarseness of the approach here, necessitated 
by the large basis sets, may result in some important levels 
not being detected. Additionally, for the remaining resonances 
with short-range character, very few possess strong coupling 
strengths to the metastable manifold. We do note that there are 



some particular resonances which stand out in that their short- 
range spin-stretch character is high with P str > 80%. It is these 
levels that we believe will be the most likely to be observed 
in experiment. We also note that the majority of resonances 
appear to be dominated both by T = 1 and by a projection of 
= 1. 

4 Conclusions 

The bound states of the fermionic 3 He(2 3 Si)+ 3 He(2 3 Pj) sys- 
tem, where j = 0,1,2, have been investigated using the re- 
cently available ab initio short-range > 3,5 £+ and 1 ' 3,5 II g ,„ po- 
tentials computed by Deguilhem et al. (8 |. Single-channel and 
multichannel calculations have been undertaken in order to in- 
vestigate the effects of Coriolis and non-adiabatic couplings. 
In contrast to the situation for the 4 He* system J9J where the 
effect of these couplings on the large number of bound lev- 
els below the lowest asymptote ( j = 2) could be studied, most 
of the levels for the 3 He* lie above the lowest asymptote and 
become resonances due to couplings to open channels. 

The single-channel long-range levels obtained in the present 
investigation differ significantly from those found by Dickin- 
son ATI , both in their patterns and energies. Dickinson reports 
nine levels for the + symmetry, 16 for 0~, six for 3 + and four 
for 3~ whereas we find 22 levels for + , 35 for 0~, seven for 
3 + and five for 3~. We also find numerous levels for the l* 
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Table 3 Energies, in units of MHz, of resonances in 3 He(2 3 Si)+ 3 He(2 3 Py) that are most likely to be observable in experiment. Energies 
given are relative to the specified asymptotic energy E^. The predissociation width r pre , short-range spin-stretched character P str , coupling 
strength is/ str and largest contributing basis of <j> are listed for each level. 



Pt 

= 1780.85 MHz 

-1 
-1 
-1 
-1 
-1 

: 6292.91 MHz 

-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 

: 6739.70 MHz 
-1 
-1 
-1 
-1 

: 8520.55 MHz 
-1 
-1 
-1 
-1 
-1 

: 13032.61 MHz 

-1 
-1 
-1 
-1 



E (MHz) 

-1283.40 

-705.27 

-301.47 

-110.16 

-71.68 

-1951.86 

-1808.40 

-1179.49 

-958.79 

-848.05 

-812.01 

-779.20 

-601.37 

-499.02 

-324.32 

-313.77 

-193.32 
-186.29 
-38.44 
-11.52 

-1029.50 

-840.44 

-513.78 

-380.48 

-245.71 

-1996.78 
-680.18 
-552.73 
-508.69 



r pre (MHz) 

15.32 
71.56 
33.13 
19.02 
5.15 

69.9 
60.4 
59.8 
60.9 
46.9 
44.2 
44.9 
55.7 
86.6 
54.4 
58.0 

58.1 
36.0 
64.7 
76.6 

20.1 
24.9 
42.3 
53.7 
55.3 

91.1 
86.2 
76.3 
92.0 



gtt (%) 

49.6 
82.0 
90.2 
44.5 
64.4 

47.3 
49.1 
51.8 
76.1 
70.8 
62.2 
76.9 
49.1 
58.0 
68.1 
75.2 

51.3 
55.1 
59.7 
58.0 

76.5 
83.8 
84.8 
84.0 
73.0 

76.6 
80.5 
74.2 
76.1 



(E h ) 



0.372 
0.177 
0.267 
0.375 
0.280 

0.278 
0.313 
0.259 
0.234 
0.234 
0.268 
0.320 
0.285 
0.337 
0.281 
0.312 

0.301 
0.372 
0.290 
0.269 

0.204 
0.212 
0.204 
0.190 
0.213 

0.157 
0.182 
0.172 
0.190 
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and 2 ± symmetries for which Dickinson could not find any 
states. These differences are not unexpected as our expression 
d29l for the matrix elements of H e \ differs from that of Dick- 
inson by an overall phase factor and the phase of the A g — A„ 
term. By using the expressions given by Dickinson and with 
some modification of the values for hyperfine structure, our 
single-channel calculations were able to reproduce the results 
of Dickinson to within 5%. 

The possible experimental observability of the theoretical 
levels has been assessed using criteria based upon the short- 
range character of each level and their coupling to metastable 
ground states. Although the bound states below the lowest 
asymptote and most of the large number of resonances above 
this asymptote do not satisfy our observability criteria we are 
able to identify some 30 resonances which are promising can- 
didates to be observed in experiment. Unfortunately, the levels 
that were found in the single-channel calculations were not 
able to be linked to any of the predicted multichannel reso- 
nances. This is because we only have information regarding 
resonances that have small predissociation rates, instead of for 
the complete set of states. Hence it is very difficult to ob- 
serve the change of behaviour of a single-channel bound level 
after the non-adiabatic and Coriolis terms are included. In 
contrast, the 4 He* calculation focused on multichannel bound 
levels which allowed a comparsion between the complete set 
of single-channel and multichannel levels. For the short-range 
levels, this lack of connection implies that the non-adiabatic 
and Coriolis couplings modify the character of the levels such 
that they are no longer observable in experiment. However, 
because the resonance search is costly to perform, we can- 
not make the same statement for the purely long-range single- 
channel levels. Hence, we also recommend that future exper- 
iments also search for the levels that are marked in Tables [TJ 
and|2]as observable. 



Appendix: Basis states and matrix elements 

The unsymmetrized body-fixed (molecular) states in the cou- 
pling scheme ( fTTT ) are 

I (7i h hf\ )a , {Yihhf2)Bj, %, T, m T ) = 

\T,m T ,a f ) LILLE 

a fl a f2 a h a h n,^ n Ll a Ll n Sl &s 2 

„rf\hf fhhh fhhh p^\S\h 
x a h a h n f ^n }l n,, a h a h n, 2 a h L n L[ n S] a h 



x|y 2 £2 L2 £Xs 2 ) B |/ 2 a i2 ) B 



(50) 



where the transformation between the molecular and space- 
fixed states is, for example, 



The states of the dimer system must be constructed to cor- 
rectly include the symmetries present in the system. Impor- 
tantly, they must be eigenstates of the total parity and nuclear 
permutation. The total parity operator Pj is equivalent to the 
action of P/P^P/A^v where the action of the inversion operators 
Pl,Ps and Pj on the orbital, electronic spin and nuclear spin 
space-fixed states respectively is 



P L \Lim Li ) A = Pi\Lim Lj ) B , 
Ps\Sim Sj ) A = \Sim Sj ) B , 
PiVm^A = \Umu)B 



(52) 



then 



and 



where P, is the parity of the atomic state. The nuclear permu- 
tation operator Xn reverses the molecular axis which is equiv- 
alent to A f>B and (0,<p) — >• (it — 9,(p + it). Noting that 

= (-l)%,-a ; (9.fl,0) (53) 
P T \T,m T ,n f ) = (-l) T \T,m T -n f ) (54) 

Pr\(a x ) A ,(ai)Bj&fJ,m T ) = 

A^2(-l) / " r |(ai)A,(a2) B ,/,-%,r,mr> (55) 

where we have introduced the notation a,- = {Yi>ji>ii>fi}- The 
eigenstates of Pj are therefore given by (IT4l . 
Since a# is equivalent to PjPiPsPi where 

P L P s Pi\(a 1 ) A ,(a 2 )BJAf' T ^ m T) = 

(-l)/i+/2-/p 1 P 2 |(a 2 ) A ,(a 1 ) S) /,a / ,r,m r ) (56) 

then the action of Xjy on the states (TT41 is 

X N \(a l ) A ,(oiQ,)B,f,^,T,mT;PT) = 
p T p x p 2 (-l) f[ +h-f+ N ^ 

x\{cb) A ,{<u{)Bj,W,m\Pr) (57) 

so that the eigenstates of Xn are dTBl . 

The relationship (fJTJi between the bases (fTTT i and (fT9b is 
obtained by first using 



l7ii2 Ll n Sl ) A |^n L2 n S2 ) B = 

L n il a z . 2 a L L a 5l a 52 n 5 



Lsa L a s 
x\( Yl ) A (n)B,LSCl L £l s ), 



(58) 



|7fly)=L^(9,0,O)Umy) 



(51) 



and expressing sums over Clebsch-Gordan coefficients as 9 — 
j symbols to give 

\(a l ) A ,(a 2 ) B J^J,m T ;P T ) = 
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\T, mT j) EEL 



f h h 






i 2 




I h h 


111 






;} 


{ f\ h 




1 /i 


72 





x|(7i)A(72)B,"£2 L £l s )|/ 1 n,- 1 ) A |/ 2 a 2 )B. (59) 
Introducing the coupling coefficients, for example, 

llsliZ = [(2L+l)(2S+l)(2j 1 + l)(2j 2 +l)]a 

U L 2 L 



h h i 



(60) 



and using, from ( |T9l , 

I(7i)a(72)b,"£2 l £2 s ) =N w (\g) + \u)) (61) 

then 

\(a 1 )A,(a 2 )B,f,<l>,T,mT;PT) = 

\T,m T ,t) E E E Cjftm, 

^%*^^(k) + l"))l»'i«ii>A|feflfa)j». (62) 

The state with A «-» B is obtained by reordering the angular 
momenta subscripted with 1 and 2 in all Clebsch-Gordan and 
9 — j symbols and using 

|(72)A(7i)B,"n L i2 s >=M VJ PiP 2 (-l) Ll+i2 - L+Sl+52 - 5 (|^)-|«)) 

(63) 

to give 

\{a 2 )A,{a l ) B ,f,<p,T,m T ;P T ) = 

\T,m T A) E E E (-l) /l+/2+/+2( - 



x(\g)-\u))\i 2 Cl i2 ) A \hCl h ) B . 



(64) 



Forming the combination ( [TBI then yields (fJTJi. 

The matrix elements of //el are diagonal in 0. Using the 
explicit states (f25T > then 



(ai,04,f |tf e i|O£i,OC2,0) = 

(T',m' T ,<l>'\T,mT,<l>) L L L L i' 1 )^' 2 
t s'n' L Q' s sn L n s 

Aij'AfiS'i'iijfxhSjz} 1 ' 2 

^ r fi'f r jif r is'j' r lSj 

X ^a' j Q.' i (j>'^SljQ, i (l>^n.' L Q.' s Q.' j ^Q. L Q, s Q, j 



1 h f ) { 1 h j 

x{ 1/2 1/2 i' W 1/2 1/2 i 

A A f J I /1 /2 / 
X \ 1 / S> J \ 1 7 5 

x^[(4(g'l+<r(«')|Hei|(4l^)+^l"))] 

x<(ri) A ,(ii) J) j / ,n;.|(i 1 ) A ,(j 1 ) B ,j ) n i ) 



(65) 



where cf r = 1 + (-l)'Pr and c|' r = 1 - (-1)'Pt. The orthog- 
onality of states reduces this to 

(«(,«2^'l^el|«l,«2,0) = 

5 ^'EE E (-l) j ' 2+ HfflfU2jfif2h} 1/2 [Si] 

vr ;"''/' r isf r isj 

A u a^. n, u n y a,- ^ a L n s n'. ^ n L % a ; 

x j 1/2 1/2 "i 1 J 1/2 1/2 i 1 { j J, j 

{ A A f J I /1 / 2 / J 1 J 

X { ! J S }^^( 2S+1 A.^W+^) 
+ cJ' r ( 2s+1 A+ + ££)] (66) 
where tj = {71 , 72, 0, T,mT,Pr}- Similarly we can show 

{a[,(4,-^>'\6 e i\(Xi,a2,-^) = 

{^\) f -f'{a[,a^A'\H s \\a Xl a 2 ^). (67) 

The phase factor (— l)-^ - ^ is cancelled by the factor arising 
from (fT8T l, thus ensuring the two contributions (l66t and (l67l 
interfere constructively. Conversion of the Clebsch-Gordan 
coefficients into 3 — j symbols finally yields ( |29l . 
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